YITP-98-12, 



y-qc/ 9904054 



Apparent Horizon Formation and Hoop Conjecture in 
Non-axisymmetric Spaces 



Takeshi Chiba* 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-01, Japan 

(February 7, 2008) 



Abstract 

We investigate the validity of Thome's hoop conjecture in non-axisymmetric 
spacetimes by examining the formation of apparent horizons numerically. If 
spaces have a discrete symmetry about one axis, we can specify the boundary 
conditions to determine an apparent horizon even in non-axisymmetric spaces. 
We implement, for the first time, the "hoop finder" in non-axisymmetric 
spaces with a discrete symmetry. We construct asymptotically flat vacuum 
solutions at a moment of time symmetry. Two cases are examined: black 
holes distributed on a ring, and black holes on a spherical surface. It turns 
out that calculating C is reduced to solving an ordinary differential equation. 
We find that even in non-axisymmetric spaces the existence or nonexistence 
of an apparent horizon is consistent with the inequality: C ^ AttM. 

PACS numbers: 04.20.Dw; 04.25.Dm; 04.70.Bw; 97.60.Lf 



Typeset using REVT^ 

*Current address: Department of Physics, University of Tokyo, Tokyo 113-0033, Japan 

1 



I. INTRODUCTION 



In 1972, Kip Thorne proposed the hoop conjecture [|I|; horizons form when and only when 
a mass M gets compacted into a region whose circumference in every direction is C ^ AnM. 
However, it took nearly a decade until the studies appeared which directly address the 
conjecture. This may be because of the vagueness of the statement, i.e., the circumference 
C of the horizon and the mass M of a body are ambiguously defined. 

For axisymmetric bodies, the definition of C is rather simple as either the minimum 
equatorial circumference or the minimum polar circumference @J^. Hence there exist a lot 
of works both by numerical and by analytical approaches addressing the hoop conjecture in 
the axisymmetric system 

For non-axisymmetric bodies, on the other hand, we didn't have any satisfactory def- 
inition of C until we proposed a definition using closed geodesies 0. In this paper we 
demonstrate how to calculate C in numerically generated spacetime without axisymmetry 
but with a discrete symmetry and assess the validity of the hoop conjecture. We study 
the formation of apparent horizons because apparent horizons can be identified locally and 
provide the practical definition of a black hole. It is important to emphasize here that event 
horizons may clothe the singularities even if apparent horizons do not appear on the spa- 
tial slices considered. The existence of an apparent horizon is a sufficient condition for the 
formation of an event horizon 0. 

Unfortunately there exists an obstacle to study the formation of apparent horizons nu- 
merically because we do not have an efficient and robust method to find apparent horizons 
in general non-axisymmetric spaces, although there has been recent progress iPHTll] . This is 
an important aspect because we are interested in the existence or nonexistence of an appar- 
ent horizon. On the other hand, if we allow spaces to have a discrete symmetry (but not 
continuous one), we can solve the equation for the surface of an apparent horizon, which is 
generally a nonlinear elliptic equation, as a boundary value problem |jlO[ . Then the numerical 
treatment of the problem becomes quite similar to that in the axisymmetric case. 

In this paper, taking advantage of recent progress both in the definition of the circumfer- 
ence and in finding apparent horizons numerically, we study, for the first time, the condition 
for the formation of apparent horizons in the light of the hoop conjecture in spaces without 
axisymmetry but with a discrete symmetry. We present the momentarily static vacuum 
configurations. We consider two families of configurations as a demonstration: black holes 
distributed on a ring, and black holes on a spherical surface. We increase the number of 
black holes keeping the total mass constant and explore the existence of a common apparent 
horizon. As for the definition of a mass we simply adopt the ADM mass, because of its 
uniqueness and definiteness at least for the initial data on a spacelike hypersurface. 

This paper is organized as follows. In Sec. 2, first we review the method for finding 
apparent horizons in spaces with a discrete symmetry, and then we also review a definition 
of the circumference of a body which has been proposed by us, and finally we present the 
results of our numerical analysis of the initial data. Section 3 is devoted to summary. 



II. INITIAL DATA ANALYSIS 
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A. Finding Apparent Horizons 



A marginally trapped surface is a closed 2-surface S where the expansion 9 of future- 
directed outgoing null vectors £^ normal to it vanishes 0. An apparent horizon is defined 
as the outer boundary of a connected component of the trapped region. When the appar- 
ent horizon forms, there always exists the event horizon enclosing it if there is no naked 
singularity and the null convergence condition is satisfied 0. 

Let be the outward-pointing spacelike unit normal to S and be the unit normal to 
a time slice. Then £^ can be written as i'^ = n'^ + and thus 

e = Dis' + Kijs's^ -K = (2.1) 

on 5*. 

In this paper, we consider the conformally flat space for simplicity 

dp = tp'^fi^dx'dx^, (2.2) 
where fij denotes the fiat metric. Then Eq. (|2.1| ) is rewritten as 



e 



^2(^2+^2)3/2 



vq cot 9 + 



W^'' ^^^^ ) (r' + rj + 

yip ip r"^ ip r"^ sin^ 6 " 




- {re^ - r^, cot Q) + -fr^^ + K,,s's' -K = 0, 

sm'' 9 ' J.2 gj^2 Q 




where r{6,(f)) parameterizes the surface of the apparent horizon and rg = dr/dO. This is 
the elliptic partial differential equation about 9 and (f). It is difficult in general to solve 
numerically Eq. (|2.3| ) efficiently and robustly, and several methods have been proposed so far 

Since we are interested in the criterion of the formation of horizons, it is important to find 
horizons robustly. In axisymmetric spaces, Eq. ( |2.3| ) can be solved as a two-point boundary 
value problem [jl2[. And it is proved to be a robust method to find apparent horizons since 



the regularity of the apparent horizon on the symmetry axis is automatically guaranteed. It 
is desirable to solve Eq. ( p.3|) in a similar manner in the non-axisymmetric space. To do so, 
we assume that the space has a discrete symmetry such that r{6, (f) + ti) = r{9, 0). We also 
assume the space has a reflection symmetry about 6 = 7t/2 plane. The latter assumption 
is just for simplicity and is not essential for the method. Then the axis 6 = becomes the 
symmetry axis as in the axisymmetric case. Thus we can specify the boundary condition at 
6 = 0,77/2 as 

re = 0. (2.4) 



Hence we can solve Eq. ( |2.3| ) as a boundary value problem |T^ although the space is not 



axisymmetric. By finite-differencing Eq.( p.3|) with taking account of the discrete symmetry 
and the boundary condition Eq.(|2.4|), it becomes a matrix equation, and we solve it iter- 



atively. The detailed numerical implementation of the method is given in |T0[. When the 
apparent horizon forms, we compute its area A 
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A = A r tp'r^ sin ed<j)deJ 1 + ^ + -^r^. (2.5) 
JO Jo V sm 6* 

In the numerical results shown below, we discretize 6 and (p as follows 

4',= (l-]:)^, (2.6) 



2 J 2Ne V 2 J N, 

with 2 = 1,..., Ai'e, and j = 1, . . . , A^",^. We typically take the grid numbers Nq = = 64 
and sometimes Nq = = 128 when a higher resolution is required. 



B. Calculating Hoop 

To gauge the hoop conjecture, we calculate the proper lengths of geodesies of various 
orientations that enclose the entire configuration. Due to the symmetry such that the re- 
flection symmetry about the equatorial plane and the discrete symmetry about the ^ = 
axis, the equatorial {6 = vr/2) plane and some =constant surfaces take the special posi- 
tion. The appropriate circumference should be the maximum of (i) the minimum equatorial 
circumference C™"^ and (ii) the minimum polar circumference along (p = constant plane 

^pol ■ 

The geodesic equation on the equatorial plane is given by 

^ (2.7) 




We note that C™™ is well-defined only if the equatorial plane symmetry is assumed as in the 
recent 3D dynamical calculations on coalescing binary NSs/BHs |jl3|. On the other hand. 



geodesies along the = constant plane do not exist in general. However, if we restrict 
ourselves to the plane with the symmetry such that = 0, then such geodesies do exist 
and follow the equation: 



_2^r| ^ / ^ + l]r^^ 2%o + (2^ + -] r\ (2.8) 



Eq.( |2.7| ) and Eq.( |2.8| ) are the same except for the replacement (p —>■ 6, and we solve them by 
the fourth-order Runge-Kutta method. We vary r and r^{rQ) at = 0(^^ = 0) and compute 
C encompassing all black holes using Eq. (PT7|) (Eq. (|2.8D) up to = ^{O = 7r/2). In this way 
we obtain the minimum circumference: 

CT,''' = 2£^'{r' + rlfd<p, (2.9) 

= 4 1^^' {r' + rl)Ue. (2.10) 

We note that when we compute the minimum polar circumference, we first minimize the 
circumference with being fixed and then maximize over which satisfies = 0. 
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C. Numerical Results 



Now we construct the vacuum initial data. We consider the time-symmetric initial data, 
that is, Kij = 0. Then the only constraint equation we have to solve is the Hamiltonian 
constraint equation 

Alp = 0, (2.11) 

where A denotes the flat Laplacian. We consider the following configurations: (1) black 
holes distributed on a ring, (2) black holes on a spherical surface. 



1. Case (1) 

We show the numerical results for the case (1): black holes distributed on a ring. 

We consider black holes (where N is an even number) of an equal mass dis- 
tributed on a ring of radius a such that each black hole's Cartesian coordinate is 
(a cos(2'7r//A^), a sin(27rZ/A^), 0), where / = 0, . . . , — 1. Then the solution of the constraint 
Eq.(pTTD is 

M ^"^ 1 
^(r,^,0) = 1 + — V ^ (2.12) 

,=0 Jr2 - 2ar sin ^ cos (0 - 2t^I/N) + a? 



where M is the ADM mass of the total system. In the limit of large A^ the configuration 
becomes the singular ring, and the solution becomes 

^{r, d) = l+ ^ . K{kI (2.13) 

where K{k) is the elliptic integral of the first kind with 

r, 4ar sin 6 

= 2.14 

r"' + 2ar sm fc' + 

The results are shown in Fig.l and Table 1. The top view of the apparent horizon as well 
as the bird's-eye view of it is shown. The shape of the apparent horizon indeed respects the 
discrete symmetry of the configuration. As A^ increases, the shape becomes pancake-like. 
The maximum of (C^*", C^oD is shown in the table. We calculate the circumference C of the 
system irrespective of the existence of an apparent horizon. Even if the apparent horizon 
forms, the appropriate circumference is generally not located on it. The circumference on 
the apparent horizon can be significantly larger [0,^,|T^ . 



When the apparent horizon forms, its area A normalized by IGvrM^ is computed be- 
cause in the time-symmetric initial data the following inequality holds if the final state is 
stationary: 

A < Aeh < Af{= 16nMf) < IGnM^. (2.15) 

Here Aeh is the area of the event horizon. We note that the existence of an apparent horizon 
implies the existence of the event horizon enclosing it 0. From this and the fact that the 
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apparent horizon is the minimum surface in the time-symmetric initial data |T^, the first 
inequahty holds. ^/ is the final stationary state black hole's area, and the second inequality 
follows from the area theorem [|^. From the uniqueness theorem [|I6|, the final state is the 
Schwarzschild black hole with mass Mf because we are considering the non-rotating system. 
Mf should not be larger than the initial mass M because gravitational radiation conveys 
the energy during the dynamical evolution. 

In the cases N = 2 and — > cxd, we also show the numerical results using a 



2D(axisymmetric) apparent horizon finder [|T^] for comparison and the check of our nu- 
merical results. We find good agreement. 

2. Case (2) 

Next we show the numerical results for the case (2): black holes on a spherical surface. 
Black holes are distributed on a sphere of radius a such that each black hole's Cartesian 
coordinate is 

{asm{Ttk/n) cos{27cl/n), asm{7ck/n) sin{2TTl/n), acos{TTk/n)), (2-16) 

where k = 1, . . . ,n — 1,1 = 0, . . . ,n — 1. We also consider two additional black holes located 
at ^ = 0, and tt on the surface so that in the large limit the configuration becomes a 
singular spherical surface. The total number of black holes is then = n{n — 1) + 2. We 
assume n is even number. The solution of the constraint Eq. is simply 

M / 1 1 \ 
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2N \ v^r^ + 2ar cos 9 + a"^ ^/r"^ — 2arcos9 + a 
2-^ k=i 1=0 Jr"^ + a? — 2ar (sin^^ sin(7rA;/n) cos{(f) — 2T[l/n) + cos^^ cos(7rA;/'n)) 



The results are shown in Fig. 2 and Table 2. For smaller N the apparent horizon is bumpy. 
Since the system becomes spherically symmetric in the large N limit, A/IQttM"^ becomes 
close to 1 as increases and C does not clearly distinguish the existence or non-existence 
of the apparent horizon. 

III. SUMMARY 

We tested the hoop conjecture in spaces without axisymmetry but with a discrete sym- 
metry. The existence or nonexistence of an apparent horizon encompassing all black holes 
is qualitatively consistent with the hoop conjecture proposed by Thorne. From Table 1 and 
Table 2, we find that if the circumference C satisfies C/4ttM ^ 1.168, then the configuration 
will be surrounded by a common apparent horizon. 

In non-axisymmetric spaces with a discrete symmetry, searching for hoops is less com- 
putationally demanding than searching for apparent horizons: just solving the ordinary 
differential equation. Further if the equatorial plane symmetry is assumed as in the recent 
3D dynamical calculations on coalescing binary NSs/BHs []T3|, C™" is then well-defined, and 



therefore even in more general spaces the hoop concept can be a useful diagnostic for the 
final fate of the collapsed object. 
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FIGURE CAPTION 



Fig.l. The shape of the apparent horizon for the case (1) near the critical separation is 
shown for (a) N ^ 2, (h) N = 4, (d) TV = 6, (d) iV = 10 and (e) = oo. The 
top view of the apparent horizon as well as the bird's-eye view of it is shown. As N 
increase the shape becomes pancake-hke. The coordinates are in units of M. 

Fig. 2. The shape of the apparent horizon for the case (2) near the critical separation is 
shown for (a) A^ = 4, (b) A^ = 14, (c) A^ = 32 and (d) A^ = 58. For smaller A^ the 
surface is bumpy. As A^ increases it becomes spherical. 



TABLE CAPTION 

Table 1. Properties of black holes distributed on a ring. Here "2D" indicates the numerical 
results using an axisymmetric apparent horizon finder. "Y" means that apparent hori- 
zon enclosing all black holes are found; "N" means the opposite. C — max(C^*", C^|"). 

Table 2. Properties of black holes on a spherical surface. 
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Table la: N=2 case 

a/M AH? C/47rM A/I&t^M'^ A/167rM^(2D) 
0.38 Y TUA (19780 (19780 
0.383 Y 1.115 0.9769 0.9769 
0.39 N 1.118 



Table lb: N=4 case 



a/M AH? C/A-kM A/ienlP 
0.41 Y TliO 0.9865 
0.42 Y 1.144 0.9841 
0.43 N 1.149 



Table Ic: N=6 case 



a/M AH? C/AnM A/WttW 
0.43 Y 054 0.9883 
0.44 Y 1.159 0.9860 
0.45 N 1.163 



Table Id: N=10 case 



a/M AH? C/AttM A/WttM'^ 
0.45 Y Ties 0.9875 
0.46 Y 1.168 0.9857 
0.47 N 1.173 



Table le: N=oo (ring) case 



a/M AH? C/AttM A/IQttM^ A/167rM^(2D) 
~0A9 Y 1.155 0.9817 0.9812 
0.495 Y 1.158 0.9794 0.9797 
0.505 N 1.163 
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Table 2a: N=4 case 



a/M AH? C/A-kM A/WnM'^ 
0.40 Y rise 0.9890 
0.41 Y 1.140 0.9870 
0.42 N 1.145 



Table 2b: N=14 case 



a/M AH? C/AttM A/IQttM^ 
0.40 Y r068 0.9982 
0.43 Y 1.076 0.9965 
0.44 Y 1.079 0.9954 
0.45 N 1.082 



Table 2c: N=32 case 



a/M AH? C/AnM A/WnTP 
0.40 Y 1.054 0.9992 
0.44 Y 1.062 0.9982 
0.45 Y 1.064 0.9975 
0.46 N 1.066 



Table 2d: N=58 case 



a/M AH? C/AnM A/lGifW 
0.40 Y r051 0.9993 
0.44 Y 1.058 0.9986 
0.45 Y 1.060 0.9983 
0.46 N 1.061 



FIGURES 
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